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If cosmological cold dark matter (CDM) consists of light enough bosonic particles that their phase- 
space density exceeds unity, they will comprise a Bose-Einstein condensate (BEC). The nature of 
this BEC-CDM as a quantum fluid may then distinguish it dynamically from the standard form of 
CDM involving a collisionless gas of non-relativistic particles that interact purely gravitationally. We 
summarize some of the dynamical properties of BEC-CDM that may lead to observable signatures 
in galactic halos and present some of the bounds on particle mass and self-interaction coupling 
strength that result from a comparison with observed galaxies. 
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I. INTRODUCTION 

Astronomical observations suggest the presence of non- 
baryonic, non-relativistic ('cold') dark matter (DM), 
comprising around 23 per cent of the energy density in 
the Universe. The particle nature of dark matter remains 
elusive, however, despite ongoing efforts to detect it di- 
rectly in search experiments or indirectly via imprints 
on astrophysical observations. In large N-body simula- 
tions of structure formation, cold dark matter (CDM) has 
been modeled as a collisionless gas, which only interacts 
gravitationally. We refer to this as standard CDM. De- 
spite the successes in reproducing the large-scale struc- 
ture, as well as reproducing flat rotation curves at the 
outskirts of galaxies, standard CDM seems to be in con- 
flict with observations of galactic small-scale properties. 
These are notably the overabundance of subhalos around 
big hosts of Milky- Way size and beyond, and the failure 
of reproducing flat cores in the centers of dark-matter 
dominated dwarf and LSB galaxies. Both features ap- 
parently contradict astronomical observations, and have 
been the subject of active research in the past decade. 
One approach has been to study the combined effects 
of the collisionless CDM component and the dissipative 
baryonic component, to see if the dissipative hydrody- 
namics of the latter can affect a cure. Another possible 
solution has been to go beyond the simplifying assump- 
tion that the DM particles are cold and/or collisionless. 
As a result, there has been a recent revival of investi- 
gations of structure formation with non-standard dark 
matter candidates, like warm dark matter e.g. Schneider 
et al.jHI, self- interacting fermionic dark matter e.g. Ahn 
& Shapiro Koda & Shapiro [l^ and Bose-Einstein- 
condensed dark matter (or scalar- field dark matter), e.g. 
Woo & Chiueh [H and Suarez & Matos [35i], all pro- 
posals of which are able to suppress the formation of 
self-bound structures below a certain scale as well as to 
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flatten central profiles, depending on the respective par- 
ticle parameters. 

In this paper, we will address Bose-Einstein-condensed 
DM (BEC-CDM). We are particularly interested in a 
class of models describing self-interacting bosonic dark 
matter, the particles of which are so light that they col- 
lectively occupy their ground state below a certain tem- 
perature, forming a Bose-Einstein condensate (BEC) in 
the early Universe. Then, this state can be well described 
by a scalar field, the so-called wavefunction of the con- 
densate. We will assume that the wavefunction respects a 
f7(l)-symmetry, such that the number of particles is con- 
served. There is thus no self-annihilation of DM in this 
scenario, in contrast to models considered for instance by 
Tkachev [s^ . Recent advances in particle theory predict 
the generic existence of bosons as light as or much lighter 
than the QCD axion - the 'classic' bosonic DM candidate, 
values ranging from about 10""^^ to > 10~^ eV/c^, which 
can serve as the CDM_in the Universe (see e.g. Giinther 
& Zhuk [Hi, CarroU [8|, Arvanitaki et al. The the- 
oretical description of these very light bosons in terms of 
scalar fields leads to halo dynamics, which can be under- 
stood as the solutions of nonlinear wave equations (e.g. 
Alcubierre et al. [2|] , Chavanis \1Q^ , Matos & Ureiia-Lopez 
(23 |. Sin [sl], Urefia-Lopez & Guzman [13]), in contrast 
to the A^-body dynamics of standard CDM. As a matter 
of fact, much less investigation has been pursued in the 
literature so far, in particular with respect to the non- 
linear stages of structure formation for this type of DM. 
Therefore, we are in need of a better understanding of 
how much of the parameter space of this form of dark 
matter is able to reproduce the successes, while resolving 
the failures, of standard CDM. 

In what follows, we will summarize some of the recent 
work in this field, highlighting some of our own, with 
apologies for the fact that length limitations prevent us 
from attempting a more comprehensive review. In Sec- 
tion |TT1 we will describe the basic equations which govern 
BEC-CDM dynamics, showing how its quantum nature 
leads to fluid behavior. We will then distinguish two 
regimes according to the strength of the particle self- 
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interaction, and state how the characteristic size and 
mass of structures that form gravitationally from this 
form of DM are related to the particle mass m and self- 
interaction coupling strength g. We also present the virial 
theorem for isolated BEC-CDM halos. In Section Hill we 
summarize how the equilibrium structure of BEC-CDM 
halos can distinguish them from halos in standard CDM, 
for halos with and without rotation. We show that halos 
in the BEC-CDM model typically rotate with enough an- 
gular momentum that, if they are of the type supported 
against gravitational collapse by self-interaction pressure 
and rotation, quantum vortices are likely to form, which 
can affect halo density profiles. We also summarize a few 
examples in which BEC-CDM can be distinguished by its 
effect on baryonic structures. In Section llVi we summa- 
rize the bounds on particle mass and coupling strength 
that follow from the requirement that halos or halo cores 
not exceed some characteristic size, if halos or their cores 
are treated as isolated, equilibrium structures. Finally, 
in Section |Vl we argue that, for BEC-CDM to reproduce 
the full range of structure scales found in our Universe 
(as in standard CDM) , we must go beyond the modeling 
of individual objects in static equilibria to account for 
continuous mass infall and its dynamical consequences. 



II. QUANTUM-COHERENT DARK MATTER 
UNDER NEWTONIAN GRAVITY 

A. Fundamental properties 

The dynamic description of dark matter within galac- 
tic halos usually involves small gravitational and veloc- 
ity fields. With regard to BEC-CDM, it has proved ad- 
vantageous to consider the non-relativistic Schrodinger- 
Poisson, or Gross-Pitaevskii-Poisson (GPP) system of 
equations of motion for the dark matter BEG wavefunc- 
tion ip, as follows 

zft-^ -— AV'-hm$V + <?IV'lV, (1) 



A$ = 47rGm|V'p. (2) 

The terms on the rhs in ([T]), which govern the evolu- 
tion, are due to the quantum-kinetic energy, the gravi- 
tational potential i", and the self-interaction of identical 
bosons. The latter has been included in the usual way in 
terms of an effective interaction potential g|V'|'*/2 with 
coupling constant (or self-interaction strength) g. The 
possibly complicated particle interactions are simplified 
this way in the low-energy limit of a dilute gas: disre- 
garding higher than 2-body interactions, the cross section 
for elastic scattering of indistinguishable bosons becomes 
constant in the low-energy limit, 

a = Snal (3) 



with the s-wave scattering length a^. The coupling con- 
stant of the effective interaction is then given by 

g = 4:TTh^as/ni. (4) 

We shall note that the above CP equation is strictly 
valid only for dilute systems, which means that as must 
be much smaller than the mean interparticle distance, 
i.e. as <C n"^/"^. Also, we are restricted to g > 
{as > 0), because models described by ([U with negative 
self-interaction coupling g are not able to provide stable 
structures, since the resulting negative pressure works in 
favor of gravity. If we assume that all of the DM is in the 
condensed state, the number and mass density of DM in 
the halo is given by n(r) = |'0P(r) and p{r) ~ mn{r), 
respectively. 

The equations can be written in fluid-like form by in- 
serting into ([l} the polar decomposition of the wavefunc- 
tion (applied early by Madelung to the free-particle 
Schrodinger equation), 

V'(r,t) = |V|(r,t)e^^("^'*) = J B^^e'S(r,t) ^ 

V m 

resulting in the momentum and continuity equations. 



p_ + p(v.V)v = -pVQ-pV$- VPs,, (6) 

|^ + V.(pv)^0, (7) 

with the bulk velocity v = hVS/m. The gradient of 

g = -h^^^p/{2m^^p) (8) 



gives rise to what is often called 'quantum pressure', an 
additional force on the rhs of equ.®, which basically 
stems from the quantum-mechanical uncertainty princi- 
ple. The particle self-interaction, on the other hand, gives 
rise to a pressure of polytropic form 

Psi = .9pV(2m2). (9) 

The quantum-kinetic term provides an important charac- 
teristic length scale, as follows: with A having a dimen- 
sion of and changing to the momentum represen- 
tation, one can easily see that the characteristic length 
is essentially nothing but the de Broglie length of the 
bosons 

L ^ XdeB ^ h/p ^ h/{mv). (10) 

Since we will consider virialized, isolated objects, it 
makes sense to use the corresponding virial velocity in 
this expression for XdeB , assuming that the particles stay 
in their Bose- Einstein-condensed state after virialization, 
which is the case for the model parameters we are going to 
encounter in this paper (see also Section |IV|. BEC-CDM 
without self-interaction, Psi = 0, was termed 'Fuzzy 
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TABLE I. Lower bound on the boson mass, provided by the 
'Fuzzy Dark Matter' regime, or TYPE I BEC-CDM, for dif- 
ferent cosmological structures 





halo mass 


size 


boson 


mass 




[Mol 


[kpc] 


[eV] 






Milky Way (MW) 


10^^ 


100 


1.066 ■ 


10" 


-25 


Dwarf Galaxy (DG) 


10^° 


10 


3.371 • 


10" 


-24 


Dwarf Spheroidal (dSph) 


10« 


1 


1.066 • 


10" 


-22 


Minihalo (MH) 


10"^ 


0.1 


3.371 • 


10" 


-21 



Dark Matter' by Hu, Barkana & Gruzinov fT6|. In this 
review, we will call it BEC-CDM of TYPE I. In this 
regime, it is the quantum-kinetic term via (jlOp which de- 
termines the equilibrium size of self-gravitating objects - 
in order for XdeB to be of a certain galactic length scale, 
the particle mass m must be sufficiently large. Table |T] 
contains different typical halo sizes and the correspond- 
ing lower limits on m. However, it has been observed 
in the previous literature, e.g . Colpi, Shapiro & Wasser- 
man |12| and Lee & Lim |2l| . and we confirm it as well, 
see eau. p^ and Section Hill that a larger mass m than 
inferred from (jlOp can result in stable structures of the 
same given size, if self-interaction is included. In fact, 
considering the case in which the last term in ([T} supports 
the system against gravitational collapse, while Q — 0, 
we arrive at the opposite regime to TYPE I, which we call 
the Thomas-Fermi regime of BEC-CDM, or TYPE II for 
short. We note that this regime goes under many names; 
Goodm an Il3| calls it 'repulsive dark matter' (RDM), 
Peebles speaks of 'fluid dark matter'. We also note 
that both regimes have been already considered for re- 
lated models by Khlopov, Malomed & Zeldovich ^13], in 
studying gravitational instabilities of a primordially pro- 
duced scalar field. 

In the following, we will make use of convenient units, 
defined as in Rindler-Daller & Shapiro (RDS) [28i] : 



h 



ruH 



= 1.066- 10" 



2h 



(RM) 



-1/2 ^ 



1 kpc 



IQS Mq 



eV, (11) 



and 



gn = r/(2pi?^) = 2iTtfR/{-iM) 



2.252 • 10 



-62 



R 



M 



1 kpcy V M 



eV cm^, (12) 



with c — 1. We described the meaning and significance 
of those parameters at length in RDS [28|. It shall be 
sufficient to re-iterate here that rrijj is the characteris- 
tic mass of a non-interacting particle whose de Broglie 
wavelength is comparable to the size of a given halo, see 



Table m It is thus the smallest particle mass possible in 
order for quantum pressure to be solely responsible for 
holding that halo up against gravitational collapse. On 
the other hand, if there are density variations in the BEG 
fiuid of scale length i?, then is the coupling strength 
for which the quantum and self-interaction pressure force 
terms are equal. 

It turns out that the TYPE II regime is a good approx- 
imation as long as g/gn S> 2 is fulfilled, as we demon- 
strated in [28j]: to determine whether a BEC-CDM halo 
of a given size R is of TYPE I or TYPE II, we have to 
compare the quantum pressure and self-interaction pres- 
sure terms in ([6]) to each other, 

I - pVQI/l - VPsi\ ~ h^/{gpR^) ~ 2gH/g « 1, 

from which the claim follows. We have shown in RDS 
(28| , equ.(46), that the radius of a spherical halo is then 
related to the de-Broglie length of the boson according 
to 



12 \gH 



1/2 



A, 



deB- 



(13) 



Hence, since g/gn > 2, i?o > \deB for TYPE II BEC- 
CDM halos. 



B. Stationary systems and virial equilibrium 

In the context of BEC-CDM in the GPP framework, 
equ.lIT])-©, stationary self-gravitating halos can be de- 
scribed by wavefunctions of the form 



^(r,^)=^.(r)e-^'^*/^ 



(14) 



where the conservation of particle number fixes fi, the 
CP chemical potential. While -0 evolves harmonically in 
time, the mass density p = m\il^s\'^ and, hence, the gravi- 
tational potential $ are time-independent. Inserting this 
^ into ([l} results in the time-independent CP equation 
with eigenvalues /i, 



2rn 



(15) 



The time- independent part ■(/'s(r) itself can be decom- 
posed as usual. 



7^,(r) = |0,|(r)e*^=W 



(16) 



with both amplitude and phase depending here on posi- 
tion only. We will omit the subscript 's' in the forthcom- 
ing analysis. Systems obeying ([15]) can be studied via 
the corresponding CP energy functional given by 



Jv 



2m 



dh. (17) 



Inserting ([16)) into ([Tfl) . the total energy can be written 
as 



E = K + W + Usi, 



(18) 
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with the total kinetic energy term 



(Vv^)Vr + / ^y'd't s Kq + (19) 



2j3„ _ 



Kq accounts for the quantum-kinetic energy and T for 
the bulk kinetic energy ol the body, which comes in the 
form of rotation, internal motion, ect. Kq has no clas- 
sic counterpart, and is absent in the classical figures of 
equilibrium studied in the previous literature. Also, Kq 
is neglected in the TYPE II regime. The other terms in 
(fT8)) are the gravitational potential energy 



W . 



and the internal energy 



SI 



2nv 



(20) 



(21) 



which is determined by the particle interactions, and 
which we have defined essentially as Usi = / PsidV with 
Psi in ®- The origin of Usi is due to the repulsive 2- 
body elastic scattering of identical bosons, equ. ([3]). The 
above energy contributions enter the scalar virial theo- 
rem of an isolated (possibly rotating) BEC halo under 
self-gravity, which reads as 



2K + W + ZUsi = 0. 



(22) 



As in classical gas dynamics, ((22)) can be derived by mul- 
tiplying the equations of motion in fluid form, equ. ([6]), 
by r and integrating the resulting equation over volumes 
which enclose the system of interest. For an isolated 
body, a derivation involving a scaling argument was pre- 
sented by Wang 



III. SIGNATURE EFFECTS OF BEC-CDM ON 
HALOS AND HALO CORES 

A. Sizes and density profiles 

The equilibrium density profiles of self-gravitating 
BEC-CDM halos are solutions of (P or dH)-® respec- 
tively, along with As an important result, they are 
universal in shape. Furthermore, BEC-CDM halos have 
a finite central density. In fact, this last feature has been 
one of the motivations in the previous literature to con- 
sider this form of DM as a solution to the cusp-core prob- 
lem of dark matter dominated dwarf galaxies (DC) and 
dwarf-spheroidal (dSph) galaxies. 

In the case of TYPE I, the density profile can only be 
determined numerically: it falls off as r^^ for large r, but 



has no compact support. The radius which includes 99 
per cent of the mass reads as 



i?99 = 9.m^/{GMm^) 



(23) 



(see Membrado, Pacheco & Sahudo [25*] for more details). 
For TYPE II, on the other hand, the equation of state 
reduces to an (n — l)-polytrope, with Psi in ^ and 
Q — Q, having the well-known spherical density profile 



— pf sinc{\/ AirGm-^ / g r), 

with sinc{x) = sin{x)/x and the central density 
corresponding halo radius is then given by 



Rq 



9 



(24) 
The 

(25) 



see e.g. Goodman jl3'| for more details. 

In both regimes, TYPE I and TYPE II, the halo profile 
and size are determined by the DM particle parameters. 
Of course, this is also true for the intermediate regime: 
Chavanis & Delfini [ll| calculate numerical solutions for 
the mass-radius relationship, R = R{M), which inter- 
polate between TYPE I and II. However, for any given 
particle model, neither (|23|) and (|25|) nor the result of 
[111 are able to re-produce the fact that R must increase 
with M as we know from astronomical observations. The 
successful fitting of galaxies data using the associated ro- 
tation curves of TYPE I and TYPE II BEC-CDM hy 
Arbey, Lesgourgues & Salati [S^] and Bohnier & Harko Q 
must thus be judged with this caveat in mind. It implies 
that it is necessary to go beyond the description of non- 
rotating halos in virial equilibrium, composed of a pure 
BEC-CDM fiuid, if this DM model is to describe galactic 
structures successfully. 



B. Rotation and shape 

In RDS [28'|, we studied the properties of BEC-CDM 
halos in the TYPE II case, once rotation is taken into ac- 
count. It is generally believed that tidal torques caused 
by large-scale structure give a halo most of its angular 
momentum in the early phases of halo collapse. This 
picture has been confirmed by cosmological N-body sim- 
ulations of the standard CDM universe, which show that 
halos form with a net angular momentum such that the 
dimensionless ratio, the so-called spin parameter. 



A 



(26) 



where L is the total angular momentum and E the to- 
tal energy of the halo, has typical values in the range of 
[0.01, 0.1] with a median value ~ 0.05 (see e.g. Barnes & 
Efstathiou jH). The degree of rotational support is thus 
very small for the CDM halos which surround galaxies. 
We will be interested in the case where the BEC na- 
ture of DM affects small-scale structure and the internal 
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dynamics of halos, while large-scale structure formation 
shall follow the ACDM model to a great extent. There- 
fore, we adopt the above range of spin parameters for 
BEC-CDM halos, too. 

We describe the effect of rotation on the structure of 
BEC-CDM halos by two approximations which are based 
upon the classic models of rotating figures of equilibrium 
(see e.g. Chandrasekhar Q). The simplest description 
assumes that the halos are Maclaurin spheroids, which 
are axisymmetric, oblate, and homogeneous. Not only is 
this model fully analytical, it also provides a convenient 
background solution to perturb in determining if and 
when quantum vortex formation is energetically favored. 
In the absence of quantum vortices, however, BEC-CDM 
is irrotational, while the Maclaurin spheroid model as- 
sumes uniform rotation. Of course, this irrotationality 
can be broken locally in the fiuid, by creating quantum 
vortices, if the amount of angular momentum exceeds a 
certain minimum, as shown below. In the limit of large 
enough angular momentum that a vortex lattice devel- 
ops, in fact, uniform rotation is a good approximation, 
and so will the Maclaurin spheroids be. More generally, 
to account for irrotationality, we also consider a second 
model, that of irrotational Riemann-S ellipsoids. Since 
the classic solution for Riemann-S ellipsoids is homoge- 
neous, however, we account for the [n — l)-polytropic 
nature of TYPE II BEC-CDM by adopting the solution 
derived by Lai, Rasio & Shapiro (LRS) [l^l for compress- 
ible Riemann-S ellipsoids, based on the 'ellipsoidal ap- 
proximation' which assumes self-similar ellipsoidal den- 
sity strata. 

We denote the semi-axes of those bodies along {x, y, z) 
as (01,02,03). Maclaurin spheroids fulfill oi = 02 > 03 
with eccentricity e = (1 — (03/01)^)^/^. Using (|22l) . we 
can determine how the (mean) radius R = (010203)^/"^ 
and the spin parameter of such a halo depend on its ec- 
centricity, see RDS 





FIG. 1. Halos rotating about the z-axis for as = 1 and 
A = 0.05: Maclaurin spheroid having e — 0.302 (left-hand- 
plot) and irrotational Riemann-S ellipsoid having (61,62) = 
(0.881; 0.797) [nght-hand-plot); see also RDS Q^. 



with i?o in (|25p and 3(61,62) a function of the eccen- 
tricities, are very cumbersome and we refer to RDS p8| . 
equ.(lOl) for more details. 

For both models, the characteristic size - the mean 
radius - depends on the particle parameters exactly the 
same way as in the non-rotating case, namely R ^jm. 
The effect of the rotation is thus only to change the over- 
all multiplicative factor of this dependence as seen from 
([271) and By fixing A = (0.01,0.05,0.1), we can 

solve A — A(e) in (^5)1 and A — A(ei,e2) in (PO)) . respec- 
tively, for the eccentricities. Fig[T] shows two illustrative 
examples of our rotating halo models. 

We were able in RDS [2^ to derive the generaliza- 
tion to eau. (p4)) for the case of the (n = l)-polytropic 
Riemann-S ellipsoid analytically, assuming the ellipsoidal 
approximation of LRS [20] • Accordingly, the ellipsoidal 
density profile reads 



pfsinc 



(31) 



R 



15 



3A3(e)(l-e2)2/3 



9 



1/2 



6 arcsin e , 

A = 1 1 

5%/5 e 



: A3(e)(l- 6^)1/2 

arcsin(e) 



1/2 



(27) 



(28) 



with the ^-parameter t = r/|iy|, a measure of rotational 
support, given by (see also LRS 20]) 



t{e) = 3/(2e2) - 1 - 3^1 - eV(2e arcsin(e)), (29) 



and ^3(e) = — 2^/1 — arcsin(e)/e'^ for this model. 
On the other hand, a compressible, irrotational Riemann- 
S ellipsoid of polytropic index n = 1 must be prolate, 
i.e. its semi-axes fulfill oi > 03 > 02, and the ec- 
centricities are given by ei = (1 — (02/01)^)^^^ and 
62 = (1 — (a3/oi)2)^/2_ jj^ ^jjg^^ case, the expressions 
for the mean radius and spin parameter 



R = Rog{ei,e2) 



-1/2 



A = A(ei,e2) 



(30) 





J, = 




A .0.01 - - 




1 = 0.05 




A =0.1 







1 2 3 4 5 6 

P ' Rchar 

FIG. 2. Density profiles of the (vortex-free) (n = 

l)-polytropic Riemann-S ellipsoidal halos having A — 
(0.01,0.05,0.1), according to equ.((3l|- The profile of the 
spherical halo, equ. (|24p . is added for comparison (solid curve). 
The densities are all normalized to pf = 1; q — qRchar = 
q[g / {4:-KGm^)Y^^ . The locus of the outer surface, where the 
density vanishes, increases with A; see also RDS [2^. 
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with 



[x^ + V^/{l-e\) + z^/{l-el)]{TT/R,)\ 



pf the central density of the elhpsoid, and Rq in (|25p . It 
is plotted for different A in FigI5] 

The virial theorem (|22|) not only provides the above 
relationships for the mean radius and spin-parameter of 
a halo, but also relates the DM particle parameters in a 
characteristic way: it can be shown that the correspond- 
ing formulae are 



1/2 



mn V8A3(e)(l-e2)2/3 
for Maclaurin spheroidal halos, and 



mn V8 V 9h 



(32) 



(33) 



for Rieniann-S ellipsoidal halos, see RDS [28|. That is, 
virialized rotating halos, according to either model, will 
lie on a straight line in (logm, log(7)-space, whose slope 
is completely determined by the eccentricities or mean 
radius of a given halo. These relationships can be found 
in Fig|3]and FiglSj respectively. From those and the con- 
dition g/gH ^ 2, it follows that, in the TYPE II regime, 
m/mH ^ 1, as it must be for XdeB ^ ^o- 

In [2^, we also determined the conditions for the for- 
mation of a central quantum vortex in a given rotating 
halo. The minimum amount of angular momentum nec- 
essary to form a singly-quantized, axisymmetric vortex 
in the center of the halo is given by 



= Nh, 



(34) 



which, for a given A and L becomes a minimum condition 
on the particle mass. The relationship is linear for both 
halo models, i.e. 



m/mH = I{ei,e2)L/LQM, 



(35) 




FIG. 3. Density profile of a spherical halo with vortex in the 
center (solid line) which transitioned from a Riemann-S ellip- 
soidal halo having A = 0.05. The profile of the unperturbed 
halo (dashed line) is added for comparison; Rchar as in Fig. 
H See also RDS [2^. 



with /(ei,e2) denoting here another function depend- 
ing only on the eccentricities (see equ.(108) and (109) 
in [28]). In order to find a criterion which is, not only 
necessary, but also sufficient, we determined the DM 
particle parameters for which the central vortex is en- 
ergetically favored, i.e. the conditions which make the 
halo with vortex have less energy than a rotating, but 
otherwise vortex- free halo. For this purpose, we con- 
sider two limiting cases: Model A describes a halo with 
a high enough angular momentum that L ^ Lqm, i.e. 
the central vortex is essentially considered to be only a 
small perturbation of the total angular momentum L of 
the halo. For this case, we use Maclaurin spheroids as 
our halo model. Model B, on the other hand, assumes a 
given halo has just enough angular momentum to form 
one quantum vortex, i.e. L = Lqm- We use the irrota- 
tional, (n = l)-polytropic Riemann-S ellipsoids for this 
case. In this case, once the central vortex is energetically 
favored, its formation takes up all of the angular momen- 
tum, leaving a spherical halo with vortex, see FigO for 
an example. 

Despite differences in the assumed halo equilibrium 
models, the conclusions are the same for Model A and 
Model B: vortex formation requires a minimum particle 
mass m > merit and a minimum particle self-interaction 
coupling strength g > gcrit- For any {m,g) pair which 
satisfies m > merit and g > gcrit at a given A, this pair 
will also favor vortex formation for any larger value of 
A. Our main results are summarized in Fig|3] and [Sj 
which are independent of halo size, respectively, while 
Table HIl and [TTTl list physical units for merit, gcrit and 
the corresponding vortex size, for the cases of a typi- 
cal dwarf galaxy (DG) and a dwarf spheroidal galaxy 
(dSph), according to TablelH Since we fix A and L — Lqm 
for the case of Model B, the respective particle param- 
eters {{m/mH)crit,{glgH)crit) are uniquely determined 
via eau. (|55|) and respectively (see also Fig. [S] for 

those numbers). It turns out that for all A considered, 
vortex formation is favored for those parameters. We 
interpret this to mean that, for these same A- values, if 
L/Lqm > 1, instead, (i.e. m/mn > (m/mH)crit, ac- 
cording to vortex formation will also be favored, 
as long as g/gn > {g / gH)crit- On the other hand, for 
halos with L < Lqm, the particle parameters must sat- 
isfy m/mn < {m/mH)cr%t and g/gn < {g/gH)crit, re- 
spectively, so they will not form vortices, but can still 
be modeled as Ricmann-S ellipsoids. (This assumes, of 
course, that g/gn ^ 2, so they are still in the TYPE II 
regime.) 

We shall comment on other previous work on vortices 
in related models. Silverman & Mallet [l^ give heuris- 
tic arguments for vortex formation, but do not derive 
the critical conditions or their consequences for the DM 
particle parameters. Yu & Morgan [40| show that vortex 
lattices canprovide flat rotation curves for galaxies. Kain 
and Ling [17[ , on the other hand, find approximate solu- 
tions for the density profile of a nonrotating, spherically- 
symmetric halo with a single vortex that contains all the 
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TABLE II. Lower bounds on the boson mass and self- 
interaction coupling strength for vortex formation in TYPE 
II BEC-CDM halos with a given spin parameter A in Model A. 
^mai-values denote the upper bounds for the corresponding 
vortex core radius. 



Dwarf Galaxy 



A 


'^crit 


[eV] 


Qcrit [eV cm^] 




[kpc] 


0.01 


1.04 • 


10-^^ 


2.30 • 10-""* 


0.03 




0.05 


1.67- 


lQ-22 


5.74 • 10-®° 


0.20 




0.10 


7.33 ■ 


10^23 


1.02 ■ 10-''° 


0.47 








Dwar 


f Spheroidal Galaxy 






A 




[eV] 


gcrit [eV cm^] 


^max 


[kpc] 


0.01 


3.30 • 


1Q-.U 


2.30 • lO"*" 


3.13 ■ 


lo-" 


0.05 


5.28 ■ 


10-21 


5.74 • 10-^^ 


0.02 




0.10 


2.32 • 


10-21 


1.02 • IQ-'^® 


0.05 





TABLE III. Lower bounds on the boson mass and self- 
interaction coupling strength for vortex formation in TYPE 
II BEC-CDM halos with a given spin parameter A in Model B. 
Cmax-values denote the upper bounds for the corresponding 
vortex core radius. 

Dwarf Galaxy 



A merit [eV] Qcrit [eV cm^] Cma:r [kpc] 

aOl 1.50 ■ 10-2^ 3.59 • 10-"° 025 

0.05 3.20 ■ 10-2^ 1.53 ■ 10-^^ 1.21 

0.10 1.69-10-2^ 3.87-10-''2 2.41 



Dwarf Spheroidal Galaxy 

A merit [eV] gcrit [eV cm^] Cma^ [kpc] 

OOl 4.75 ■ 10-2^ 3.59 ■ lO"''" 002 

0.05 1.01 ■ 10-21 1.53 ■ 10-*^° 0.12 

0.10 5.34-10-22 3.87-10-^1 0.24 



angular momentum for the TYPE II case (as in Fig. [3l 
solid curve). Their estimates for the viable parameter 
space of particle mass for the Andromeda galaxy is in 
agreement with our more precise results for the Milky 
Way, see RDS [H. 

Our results strongly suggest that vortices can only 
form in BEC-CDM with a sufficiently large, positive self- 
interaction coupling strength. Vortex formation in axion 
DM without self-interaction^ which has been claimed by 
Sikivie & Yang [3l|, seems thus not viable according to 
our results. However, for TYPE II halos, i.e. those for 
which m/mH ^ 1 and g / gu ^ 2, vortex formation is fa- 
vored for the range of halo spin parameters found in the 
CDM model, for a large portion of the BEC particle pa- 
rameter space. For example, for A — 0.05, Model B yields 
{m/m.H)crit = 9.49 and {g/gH)crit = 68.00, while Model 
A yields {m/mH)crtt = 49.52 and {g/gH)crit = 2549.24 
(see Fig. |3]and[S]). Since the presence of a vortex causes 
the DM density profile to drop inside the vortex core, 
diminishing it altogether at the very center, the effect 
of vortices must be seriously considered whenever the 
TYPE II regime is studied as a model to describe galactic 
halo dynamics. At the critical values for vortex forma- 



tion and a little above, the effect of the vortex could be 
observable, as the ^maa;-values in Table HI] and HII] sug- 
gest. However, while it is true that the vortex becomes 
increasingly favored for large g (at fixed A), its size rela- 
tive to the halo size decreases, i.e. for large enough g the 
influence of the vortex will yet again diminish. 



C. Influence on baryonic substructures 

The foregoing sections exemplified some distinctive 
characteristics of Bose-Einstein-condensed DM models 
as compared to standard coUisionless CDM. A further 
interesting topic is the detailed dynamics of baryonic 
substructure within galactic halos made of BEC-CDM. 
While this may likewise provide evidence for the pres- 
ence or absence of this form of DM, it has hardly been 
considered in the previous literature. 

Goodman [l^ seems to be the first to discuss the conse- 
quences of TYPE II, as a superfluid, for a rotating galac- 
tic bar. This discussion has been continued in Appendix 
B of Slepian & Goodman [3^], with the result that the 
drag on the bar can be smaller than estimated for stan- 
dard CDM, providing a potential remedy to a problem 
which plagues the latter. 

Lora at el. [l^ , on the other hand, have recently stud- 
ied the dynamic survival of cold gas clumps and globular 
clusters in the dSph galaxies of Ursa Minor and Fornax, 
respectively, in both TYPE I and TYPE II regimes. The 
survival of those structures within the lifetime of these 
galaxies of known mass provides valuable constraints 
on the allowed boson mass and self-interaction coupling 
strength. For TYPE I, their calculations favor a boson 
mass between 0.3 - 10-^^ eV < m < 10^^^ eV, while m 
can be larger for TYPE II at fixed coupling strength, 
their constraint reading g/m^ < 7.4 - lO^^"^ eV/cm"^ (ex- 
pressed in our units), corresponding to a core size of 
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= 0.1 





No Vortex 










/'y / Vortex ^...^ 












No Vortex " 
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FIG. 4. Model A: Curves within which vortex formation is 
energetically favored in the dimensionless BEC-CDM particle 
parameter space (jn/mH , g / Qh) for halos with spin param- 
eter A = 0.01,0.05,0.1, respectively. Straight lines for the 
same A- values: halos fulfilling virial equilibrium, according 
to equ.(131|; see also RDS 
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about R > 0.7 kpc. The first result can be understood 
given the limits in Table HI while the second result is in 
accordance with previous studies, including our own ones 
as presented in Section |TT1 In both cases, Lora et al.'22'| 
find viable parts of the parameter space of BEC-CDM 
which can explain the observations, while others can be 
definitely excluded. 

The study of the dynamics of baryonic 'test bodies' 
like globular clusters or giant molecular clouds, which are 
large enough to feel the influence of the subtleties in the 
DM distribution, while still being much smaller than the 
DM halo, will provide useful constraints on BEC-CDM 
for our own Galaxy. There is valuable and interesting 
work to be done in future studies of this problem. 

Finally, we note that galaxies are often observed to 
harbor supermassive black holes at their centers, with as- 
sociated quasar luminosity which supports the idea that 
the black hole mass grew primarily by baryonic accre- 
tion. One might ask if BEC-CDM is consistent with 
this phenomenon, since its fluid behavior might imply 
a much higher accretion rate than that of coUisionless 
CDM particles (c.f. Shapiro & Teukolsky [301]). The 
gravitational collapse of non-selfinteracting scalar field 
dark matter onto a central Schwarzschild black hole in 
a spherically-symmetric space-time has been considered, 
for example, by Barranco et al. Q and references therein, 
to determine if it is possible for DM halos in this model 
to survive for cosmological time scales. As we discussed 
in Section flllBl however, halos have angular momentum. 
For accretion onto a central black hole to occur, this an- 
gular momentum must be overcome, which usually de- 
pends upon a dissipative process involving some form of 
viscosity to transfer that angular momentum outward. 
As a frictionless superfluid, BEC-CDM halos, however, 
cannot do this. Moreover, if a vortex is formed at the 
center, the lower density there will further inhibit accre- 
tion. 









(1595.07,44.58) - 


-jt — (68.00,9.49) 


1 = 0.01 


1 = 0.05 


I (17.20,5.01) 


1=0.1 - 


particle param. for which L=LQjyi * 



IV. BOUNDS ON PARTICLE MASS AND 
COUPLING STRENGTH 

We have seen that the equilibrium size of isolated, 
hydrostatic BEC-CDM halo structures is predominantly 
governed either by quantum pressure, equ.®, or by the 
poly tropic pressure due to the collective, repulsive self- 
interaction of the DM particles, equ.®, for g S 9h or 
g > gn, respectively. By opposing gravity, quantum 
and self-interaction pressures each prevent structure from 
forming on small scales, thereby imposing lower limits on 
the size of structures we can expect to find in a BEC- 
CDM universe. The lower limit set by quantum pressure 
is XdeB in (fTO]) . evaluated using the halo virial velocity 
(which also characterizes the bulk mass motions at its 
formation time). Since the smallest halos also have the 
smallest virial velocities, the most stringent limit results 
if we require XdeB ^ R for the smallest halos. This, 
in turn, imposes the lower limit on m, see Table ID For 
a given m that satisfies this lower limit, the coupling 
strength g must not exceed the value such that the radius 
of the polytrope supported by self-interaction pressure, 
Ro in or the respective generalizations for rotating 
halos in (l?fl) and ([50)1 . exceeds the size of the smallest- 
scale structure. 

The same requirement that the characteristic size Rq 
not exceed the size R of the smallest halos can also be 
used to place an upper limit on particle mass m, if we can 
establish an upper bound on the 2-body scattering cross 
section per unit particle mass, cr/m, from some other ar- 
gument. Slepian and Goodman [34] have argued that 
upper bounds on cr/m for the elastic-scattering particles 
in the SIDM model, based upon comparing that model 
to astronomical observation, should apply to BEC-CDM, 
as well. The interpretation of the Bullet cluster observa- 
tions, for example, as a nearly coUisionless merger of two 
cluster-sized halos has been found to limit cr/m for SIDM 
halos to {(j/m)raax < 1-25 cm^/g, according to Randall 
et al.d^]. We can relate a/m for TYPE II BEC-CDM to 
the characteristic size Rq in equ. (|25l) according to 



m 



8G^ 



(36) 



as pointed out by Slepian and Goodman [3J|. In order 
for (|36p to be applicable, the particle mass must be such 
that m/rriH ^ 1, in which case RDS [2^ noted that 
can be rewritten in fiducial units related to the size R 
and mass M of a given halo or halo core, as follows: 



800 1000 1200 1400 
S'Sh 



— = 2.094 • 10" 

TO 



-95 



\mH J 



FIG. 5. Model B: Dimensionless BEC-CDM particle param- 
eter space {ra/mH,g/gH) (no log-scale !) for halos with 
spin parameter A = 0.01, 0.05, 0.1, respectively, fulfilling virial 
equilibrium, according to eau. p3)) . For each A, dots denote 
those critical particle parameters for which L = Lqm; see also 
RDS (28,). 



R 



1 kpc 



3/2 



M 



108 Mq 



-5/2 2 

— , (37) 



(where 'g' here means 'grams', not coupling strength). 
This shows that a/m is very much smaller than 
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{cr/m)max ~ 1 cm^/g, unless the particle mass is 
many orders of magnitude larger than the lower bounds, 
m/niH = 1, in Table HI We can place an upper bound on 
m, in fact, if we replace Rq in eau.([551) by i?, the small- 
est halo (or halo core) size which must be produced, and 
replace a/m by {cr/'m)max, to write 

If we take as our fiducial units R — 1 kpc, and 
{a/m),nax = 1, this gives 

m < 9.193- 10"'' eV/c^. (39) 

To be self-consistent, we must check if our assump- 
tion is valid that BEC-CDM remains a pure condensate 
without thcrmalizing during virialization. If the relax- 
ation time for particle collisions to establish thermody- 
namic equilibrium at the halo virial temperature is less 
than a Hubble time (~ 10"'^'' sec), our assumption would 
break down. In that case, Slepian and Goodman sj] de- 
termined that BEC-CDM halos in the TYPE II regime 
would have cores surrounded by isothermal envelopes of 
non-condensate, which would yield density profiles in 
disagreement with observed halos (unless a/m exceeds 
(a/m) max ^ 1 cm^/g by orders of magnitude). It can 
be shown, however, that if {a/m)„iax ~ 1 cm^/g, then 
the relaxation time for achieving thermodynamic equi- 
librium is, indeed, more than a Hubble time, so thermo- 
dynamic equilibrium is not achieved for particle masses 
which obey inequality (|39|) . Hence, this upper limit is 
a self-consistent one. Apparently, there is quite a large 
range of particle mass allowed between these upper and 
lower limits. 



V. BEYOND THE POLYTROPIC SIZE LIMIT 

Nevertheless, this requirement that Rq < R for the 
smallest halos suggests there is a problem for this model 
in its simplest form, if observations require us to accom- 
modate the formation of objects as small as the smallest 
dwarf spheroidal galaxies, while at the same time serv- 
ing to explain the flattening of the density profiles in the 
cores of much larger gala xies. The remedy suggested by 
Slepian & Goodman [3J], where BEC-CDM cores in the 
TYPE II regime are enshrouded by isothermal envelopes 
of non-condensate, was unfortunately shown not to work 
by the same authors. 

We envisage a different scenario to overcome the size 
limit given by the equilibrium model as follows: we as- 
sume that ([^5]) characterizes the size of the inner core of 
larger halos, which grow larger as a result of continuous 
infall at the time of halo formation. Let us sketch this 
here in more detail. It can be shown that a classic, spher- 
ical top-hat model for the collapse and virialization of a 
cosmological density perturbation applied to BEC-CDM 
in the Einstein-de Sitter universe yields a post-collapse 



virialized object with radius Rth and uniform density 
Pq, given by 

^-^and^.fl)', (40) 

rta 3 Pta V2/ 

where rta and pta are the radius and density of the top- 
hat at the time of turn-around (i.e. maximum expan- 
sion). Then, the (virial) radius is determined from ([2^ 

as 

Not surprisingly, this Rth depends upon g and m in 
the same way the polytrope radius Rq does, so the im- 
plication is that halos of different mass must form from 
cosmological fiuctuations that collapse at different times, 
in order that pta is different. Unfortunately, more mas- 
sive halos require higher pta, since Rth is independent 
of halo mass, but higher pta requires collapse at earlier 
times, which reverses the hierarchical structure formation 
history expected for CDM. 

The crucial idea in overcoming this undesirable feature 
is the realization that eau. (|4ip neglects any internal ki- 
netic energy of the virialized object. In fact, it can be 
shown that the virial radius grows beyond Rth, once an 
(effective) kinetic term has been added, even though the 
regime of TYPE II is retained. Determining the form 
and the physical meaning of this additional kinetic part, 
and how it can advance the above description to provide 
a successful model for halo formation and structure will 
be the subject of a forthcoming paper. An immediate 
observation is the fact that the effective kinetic energy 
Kgff should not simply be a constant factor times W or 
Us I, since this will only increase the virial radius by a 
fixed factor times Rth, as can be easily shown. In fact, 
we have seen that the inclusion of uniform rotation has 
had exactly this effect, see eau. (|27l l30l) . Instead, we will 
pursue the following idea: the BEC-CDM fluid must un- 
dergo oscillations during the process of virialization due 
to its inherently quantum-mechanical nature. In fact, 
the results of work by Khlopov at al. [l8| and Guzman & 
Ureha-Lopez [l^ lend support to this idea. 

Suppose Keff describes wave motions. For the gen- 
eral argument outlined here, it is sufficient to consider a 
non- vanishing bulk velocity v as before, i.e. K^ff — T = 
J ^-v^dV > 0. Since po — const., we assume that the 
gross average in the form of T = ^Ma^ with velocity dis- 
persion cr^ = ^(v^), will capture the overall kinetic con- 
tribution due to wave motions. Additionally, in order for 
the virial radius to depend on halo mass and time of col- 
lapse tcoii, as they do for standard CDM, we require the 
top-hat density to be a fixed fraction of the background 
density at tcoii, that is, po — Apb^coih with the constant 
A depending on the background cosmology, e.g. A ~ 178 
for standard CDM in a flat, matter-dominated universe. 
The corresponding total mass is then M = ^^TR^^ Apb^coU, 
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which yields a radius 

Wicon)= . (42) 

Inserting this expression for R into (|22p with Kq — 
and solving for the unknown velocity dispersion, we get 

15 zm^ 

(43) 

To interpret this result, we observe the following; For 
(7 = 0, the last term vanishes and we recover the standard 
case, but for increasing 5 > 0, this term makes smaller. 
We can calculate the minimum mass for which (t„ = 0, 
resulting in 




= ^T^Api^^ouRTHfi (44) 



with Rth,o given by (|4T]) . Thus, the smallest mass halo 
has the minimum size of Rth,o by construction. The 
particle parameters may now be chosen such that Mmin 
corresponds to the smallest observed galaxies, as well as 
to the cores of large galaxies, e.g. Mmin — 10® Mq. 
Larger halos then follow the relationship (|42l) and have 
a non- vanishing velocity dispersion due to internal wave 
motion, according to This guarantees that BEC- 

CDM halos of mass M > Mmin would share the mass- 
radius relation of halos in the standard CDM model, if 
halos of a given mass M typically collapse at the same 
time as they do for standard CDM. 
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